The abrupt changes in the photospheric magnetic and 
Lorentz-force vectors during six major neutral-Une flares 

G.J.D. Petrie 

National Solar Observatory^ 950 N. Cherry Avenue, Tucson, AZ 85719 

ABSTRACT 

We analyze the spatial and temporal variations of the abrupt photospheric 
magnetic changes associated with six major flares using 12-minute, Ql'.b pixel"^ 
vector magnetograms from NASA's Helioseismic and Magnetic Imager instru- 
ment (HMI) on the Solar Dynamics Observatory satellite. The six major flares 
occurred near the main magnetic neutral lines of four active regions, NO A A 
11158, 11166, 11283 and 11429. During all six flares the neutral line held vec- 
tors became stronger and more horizontal, in each case almost entirely due to 
strengthening of the horizontal fleld components parallel to the neutral line. In 
all six cases the neutral line pre-flare fields were more vertical than the reference 
potential fields, and collapsed abruptly and permanently closer to potential field 
tilt angles during every fiare, implying that the relaxation of magnetic stress 
associated with non-potential tilt angles plays a major role during major flares. 
The shear angle with respect to the reference potential fleld did not show such 
a pattern, demonstrating that flare processes do not generally relieve magnetic 
stresses associated with photospheric magnetic shear. The horizontal flelds be- 
came signiflcantly and permanently more aligned with the neutral line during the 
four largest flares, suggesting that the collapsing fleld is on average more aligned 
with the neutral line than the pre-flare neutral line field. The vertical Lorentz 
force had a large, abrupt, permanent downward change during each of the fiares, 
consistent with loop collapse. The horizontal Lorentz force changes acted mostly 
parallel to the neutral line in opposite directions on each side, a signature of 
the fields contracting during the flare, pulling the two sides of the neutral line 
towards each other. The greater effect of the flares on fleld tilt than on shear 
may be explained by photospheric line-tying. 

Subject headings: magnetohydrodynamics: Sun, solar magnetic fields, solar pho- 
tosphere, fiares 
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1. Introduction 

Solar Cycle 24 has produced several major flares since NASA's Helioseismic and Mag- 
netic Imager (HMI) instrument (Schou et al. 2011) on NASA's Solar Dynamics Observatory 
(SDO) satellite (Pesnell et al. 2012) began observing continuously in March 2010. HMI 
vector magnet ogr am series covering some of these major flares have been released to the 
community since late 20lfl SDO/HMI produces full-disk vector magnetograms with 0'.'5 
pixels every 12 minutes. From flltergrams in six polarization states at six wavelengths on the 
Fe I 617.3 nm spectral line, images for the Stokes parameters, /, Q, U and V are derived, 
which are inverted for the magnetic vector components by the Very Fast Inversion of the 
Stokes Algorithm (VFISV) code (Borrero et al. 2010). The 180° azimuthal fleld ambiguity 
is resolved using the "minimum energy" method (Metcalf 1994, Leka et al. 2009). In this 
paper we analyze six 12-hour series of vector magnetograms covering six major flares occur- 
ring near the main magnetic neutral lines of four active regions, NOAA 11158, 11166, 11283 
and 11429. We describe the abrupt and permanent fleld changes that occurred during each 
flare and characterize the associated Lorentz force vector changes near the main neutral line 
of the region and within the neighboring sunspots. Here a change is deemed "permanent" if 
its effects last until at least several hours after the flare. 

Abrupt photospheric fleld changes have been observationally linked to flares in the past 
two decades; see the discussions in Sudol and Harvey (2005) and Wang (2006). Wang (2006) 
found an unshearing movement parallel to the neutral lines in flare-related longitudinal 
magnetic fleld changes in all flve 5-spot flares that he studied, implying an overall release 
of shear, but that the two polarities converged towards the neutral line during some events 
and diverged during others. Wang and Liu (2010) studied 11 X-class flares for which vector 
magnetograms were available, and found in each case an increase of transverse fleld at the 
polarity inversion line. Wang et al. (2012), Sun et al. (2012) and Petrie (2012) analyzed the 
HMI vector data for the 2011 February 15 X2.2 flare, and found similar behavior, as did 
Liu et al. (2012) for the 2011 February 13 M6.6 flare. The HMI vector data for these two 
major flares from AR 11158 have already been studied in several papers using a variety of 
methods (Wang et al. 2012, Gosain 2012, Sun et al. 2012, Liu et al. 2012, Petrie 2012, Jing et 
al. 2012). Sun et al. (2012) calculated nonlinear force-free fleld models for the coronal fleld 
from the HMI vector measurements and argued that the increase in magnetic shear observed 
at the photosphere is localized at low heights and the shear decreases above a certain height 
in the corona (see also Jing et al. 2008). Petrie (2012) found an increase in strength of the 
fleld vector at the neutral line at the time of the flare, particularly its horizontal component 
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parallel to the neutral line, accompanied by a large, abrupt, downward vertical Lorentz force 
change and a horizontal Lorentz force change acting in opposite directions on each side of 
neutral line, with the two sunspots at each end subject to abrupt torsional un- twisting forces. 
The downward and un-shearing forces were consistent with a collapse and contraction of 
fields near the neutral line. These observations support the coronal implosion interpretation 
(Hudson 2000, Hudson, Fisher and Welsch 2008, Fisher et al. 2012) where, after a coronal 
magnetic eruption, the remaining coronal field contracts downward resulting in the field 
become more horizontal at the photospheric level. Petrie and Sudol (2010) analyzed one- 
minute GONG longitudinal magnetograms covering 77 fiares of GOES class at least M5 and, 
exploring the relationship between increasing/decreasing longitudinal fields and azimuth and 
tilt angles at various positions on the disk, found that the overall distributions of longitudinal 
increases and decreases at different parts of the disk was found to be consistent with Hudson, 
Fisher and Welsch's (2008) loop-collapse scenario. Fletcher and Hudson's (2008) physical 
description of fiaring field changes remains the only detailed explanation of how a coronal 
event could cause permanent change in the photospheric field. 

The goal of this paper is to use the high-cadence HMI vector data covering six major 
fiares to extend and clarify the above results. For example, if the transverse field component 
generally increases near neutral lines during fiares, do the magnetic shear changes also follow 
a general pattern, and, if not, why not? We expect a fiare to involve the relaxation of 
magnetic stresses built up during the preceding hours and days. Can the photospheric field 
measurements shed light on how magnetic stresses are relieved during fiares? 

The paper is organized as follows. In Section [2] we will present the vector fields observed 
by HMI before and after the main fiare-related field changes took place, discussing the 
differences between these vector fields in each spatial dimension and plotting the vector 
field evolution in time. In Section |3] we will discuss the magnetic changes with reference to 
potential fields. Section [4] will describe the associated electric current changes that occurred 
during the fiares. We will derive the accompanying Lorentz force changes in Section [5l We 
will conclude in Section [HI 

2. The magnetic field vector changes 

The six fiares and their published GOES times and classes and NOAA active region 
numbers are given in Table [T] The vector field measurements covering these fiares were 
released by the HMI team in the form of 12-minute vector magnetogram images (5^, Bq^ B^) 
in heliographic coordinates {r^O^cj)) on grids with pixel size 0.03°. Figure [l] shows spatial 
maps of the vertical magnetic field component, 5^, before the fiares with the correspond- 
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Fig. 1. — The vector magnetic field before each of the six fiares. The vertical field compo- 
nent, 5^, is indicated by the color scale and the horizontal component by the arrows, with 
saturation values ±1000 G. Red/blue coloring represents positive/negative vertical field. The 
black rectangles mark the regions of major field change near the neutral lines that are used 
in subsequent analysis. The solid and dotted contours indicate strong > 1000 G) and 
quite strong {\Br\ > 100 G) fields, respectively. 
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Fig. 2. — The vector magnetic field changes during each of the six fiares. The vertical field 
change component, SB^^ is indicated by the color scale and the horizontal component by 
the arrows, with saturation values ±300 G. Red/blue coloring represents positive/negative 
vertical field changes. The black rectangles mark the regions of major field change near the 
neutral lines that are used in subsequent analysis. The solid and dotted contours indicate 
strong {\Br\ > 1000 G) and quite strong {\Br\ > 100 G) fields, respectively. 
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ing horizontal field, = {B^i^^—Bq)^ indicated by arrows. In all cases the active region 
field structure is complex but includes a prominent, highly sheared magnetic neutral line, 
indicated in each plot by a black rectangle. These rectangles were chosen balancing the 
following priorities: to align the each rectangle satisfactorily with an approximately straight 
section of each neutral line shown in Figure [T| and to capture and characterize the signifi- 
cant horizontal and vertical field changes near the neutral lines shown in Figure [2] Different 
rectangles were tried and the results were stable, although some of the patterns in the results 
became stronger or weaker for different choices of rectangle. Regular rectangular domains 
were chosen because they define directions parallel and perpendicular to the neutral lines 
and because the estimates of Lorentz force change in Section [5] are expected only to apply 
to well-resolved field changes in major magnetic structures (Petrie 2012). 

It is in this central portion of the active regions near the main magnetic neutral lines 
that most of the organized magnetic changes occurred during the fiares, as the difference 
maps in Figure [2] show. If we have observations of the photospheric vector field at two 
times, t = before the field changes begin, and t = 6t after the main field changes have 
occurred, the magnetic vector changes due to the fiare can be represented by the difference 
SB = B^St) — B(0). Figure [2] shows spatial maps of the vertical magnetic field change, SB^^ 
with the horizontal field changes (^B^ = {SB^j^^ —SBq) indicated by arrows. Each difference 
map was constructed by subtracting the last pre-fiare 12-minute image from the image with 
time stamp 24 minutes later. This is analogous to the differencing of 10-minute averages 
of MDI one-minute data by Sudol and Harvey (2005), which these authors used to validate 
their GONG field change results. These authors showed that most field changes occur during 
the first 10 minutes or so of major fiares, thus most of the field changes are accounted for in 
the difference calculation. The photospheric plasma acts continuously on the field, while the 
coronal field is believed to evolve through series of nearly force-free equilibria, punctuated 
by brief dynamical episodes such as the major fiares studied here. The major permanent 
photospheric field changes that occur during fiares correspond to the abrupt and permanent 
restructuring of the coronal field due to the fiare. 



Table 1: Flares studied in this paper. 



Date (UT) 


GOES start 
Time (UT) 


GOES Peak 
Time (UT) 


GOES End 
Time (UT) 


GOES 
Class 


NOAA 
Number 


Location 
on Disk 


2011 Feb 13 


1728 


1738 


1747 


M6.6 


11158 


S20E05 


2011 Feb 15 


0144 


0156 


0206 


X2.2 


11158 


S20W10 


2011 Mar 9 


2313 


2323 


2329 


XI. 5 


11166 


N08W11 


2011 Sep 6 


2212 


2220 


2224 


X2.1 


11283 


N13W18 


2011 Sep 7 


2232 


2238 


2244 


XI. 8 


11283 


N14W31 


2012 Mar 7 


0002 


0024 


0040 


X5.4 


11429 


N18E31 
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The difference maps in Figure [2] for the 2011 February 13 M6.6, 2011 September 6 
X2.1 and 2011 September 7 XI. 8 flares show that the vertical changes were mostly posi- 
tive/negative on the posit ive/negative side of the neutral line, weakening the vertical fleld 
on both sides of the neutral line. The map for the 2011 February 15 X2.2 flare shows the 
opposite pattern and the map for the 2011 March 9 and 2012 March 7 X5.4 flares show a 
less organized mixtures of positive and negative changes in the vertical fleld. The vertical 
fleld changes therefore do not show a pattern that generalizes across the six-flare data set. 
In contrast, the horizontal changes point in approximately the same direction as the fleld 
itself near the neutral lines of all of the six flaring region, strengthening the horizontal fleld 
there. Spatial maps of the changes in the fleld tilt angle 7 = tan~-^([5^ + and 
the total fleld strength B = [B^ + B^ + B'^y^'^ during the flare (not shown) indicate clear 
increase in the tilt angle and strength of the vector fleld near the neutral lines during the 
flares. We will discuss these increases in fleld strength in more detail below. 

As the diflFerence maps in Figure [2] show, solar flares are extremely complex physical 
events and can involve diverse magnetic changes occurring throughout the flaring region, 
most notably near major neutral lines and in sunspots. Here our goal is to shed light 
on the changes occurring near neutral lines. We now discuss the temporal proflles of the 
neutral line magnetic changes, shown in Figures [3]j6) These and subsequent plots of temporal 
changes were derived by calculating area integrals of the fleld components over the chosen 
photospheric areas in the 60-image 12-hour series, represented by the black rectangles in 
Figures [T] and [2] Figure |3] shows the evolution of the integrated magnetic fleld strength 

= / B dA, (1) 
near the neutral line over each of the 12-hour time intervals. The area A^i^ corresponds 

Table 2: Summary of changes in magnetic tilt and shear angles. The tilt angle is calculated 
with respect to the local vertical direction and the azimuthal angular displacement with 
respect to the direction normal to the neutral line. The shear and dip angles are calculated 
with respect to the reference potential fleld. 





Tilt Angle 


Azim. Angle 


Shear Angle 


Horiz. Shear 


Dip Angle 


Date (UT) 


Inc/Dec 


Inc/Dec 


Inc/Dec 


Inc/Dec 


Inc/Dec 


2011 Feb 13 


Inc 








Inc 


2011 Feb 15 


Inc 


Inc 


Inc 


Inc 


Inc 


2011 Mar 9 


Inc 


Inc 






Inc 


2011 Sep 6 


Inc 


Inc 






Inc 


2011 Sep 7 




Inc 






Inc 


2012 Mar 7 


Inc 


Inc 


Inc 


Inc 


Inc 
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Fig. 3. — Shown here are the integrated magnetic field strength B^^ (sohd fines) and B^^ 
(dashed fines) near eacfi neutral line plotted against time. Tfie areas of integration are 
indicated by the black rectangles in Figure [T] The vertical lines represent the GOES flare 
start, peak and end times. 
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Fig. 4. — Shown here are the total vertical magnetic flux B^^ near each neutral line plotted 
against time. The red/blue solid/dashed lines represent positive/negative flux. The areas 
of integration are indicated by the black rectangles in Figure [T] The vertical lines represent 
the GOES flare start, peak and end times. 
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Fig. 5. — Shown here are the integrated horizontal field parallel to each neutral line, B'^^^ 
near each neutral line plotted against time. The red/blue solid/dashed lines represent posi- 
tive/negative field, i.e., approximately westward/eastward field. The areas of integration are 
indicated by the black rectangles in Figure [T] The vertical lines represent the GOES fiare 
start, peak and end times. 
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Fig. 6. — Shown here are the integrated horizontal field perpendicular to each neutral line, 
S^^, near each neutral line plotted against time. The red/blue solid/dashed lines represent 
positive/negative field, i.e., approximately northward/southward field. The areas of integra- 
tion are indicated by the black rectangles in Figure [Tj The vertical lines represent the GOES 
fiare start, peak and end times. 
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to the rectangular region near the neutral line of each flaring region in Figures [T] and [2j 
Figure 3 also shows the equivalent integral B^^ of 5^, where 5^ = {B^ + B'^y^'^. Near each 
neutral line, according to Figure [sj the average fleld strength B^^ increased abruptly and 
permanently during each flare because of an increase in the horizontal fleld component B^^ 
there, also shown in Figure [3| where the two quantities are seen to track each other closely 
in every case. 

The behavior of the integrated vertical flux B^^ was more diverse in character. Accord- 
ing to Figure [4| around the times of the 2011 September 6 and 7 X-class flares the vertical 
flux increased somewhat abruptly near the neutral line of AR 11283. In all other cases the 
vertical fleld changes were either small or gradual. In all cases the vertical changes were less 
signiflcant than the horizontal changes. From the temporal plots it is therefore clear that 
the changes in the horizontal fleld dominated the behavior of the total fleld strength near 
the neutral line around the time of each flare, strengthening it. However, the vertical fleld 
changes did have a bearing on the physics of the 2011 September 7 XI. 8 flare, as we will see 
in Section [3l 

Figures [5] and |6] show the evolution of the integrated horizontal magnetic fleld compo- 
nents parallel and perpendicular to the neutral line, B^^ and B^^. These directions are 
deflned by the black rectangles in Figure [ij The parallel direction is the direction of the 
long edges of the rectangles, pointing approximately west. The perpendicular direction is 
the direction of the short edges of the rectangle, pointing approximately north. The flelds 
pointing in the positive (north or west) and negative (south or east) directions are plotted 
separately. The parallel flelds had a clearly dominant direction in each case and this fleld in- 
creased during all six flares. The perpendicular fleld changes were generally not as signiflcant 
as the parallel changes: contrast Figures [5] and (G) The perpendicular fleld strength increased 
abruptly and permanently during the 2011 February 15 X2.2 flare but even in this case the 
perpendicular increase was dwarfed by the parallel fleld increase. The pre- and post-flare 
evolution of the horizontal fleld was more steady in the parallel than in the perpendicular 
component and the flare-related changes much more signiflcant. 

Flare-induced line proflle changes can produce signatures that do not indicate real 
changes in the magnetic fleld as discussed by Sudol & Harvey (2005) - see the bottom 
part of their Figure 1. Working with GONG 1-minute longitudinal fleld images, Sudol & 
Harvey (2005) and Petrie & Sudol (2010) identifled artifacts by fltting a tan~^ step-like 
function to the time proflle of each pixel, applying selection criteria based on the quality of 
the function flts, and inspecting the results for representative pixels by eye. For the HMI 
vector data this approach is not as useful because the HMI vector fleld inversions are not as 
sensitive as the GONG data. 
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Because the difference images in Figure [2] are derived from measurements taken before 
and after the main field changes have taken place, thereby excluding the times when most of 
the fiare artifacts are expected to have occurred, these diflFerence images are not expected to 
be significantly compromised by fiare emission artifacts. The time series plots in Figures |3}j6] 
do not avoid these times when most of the fiare artifacts are expected to have occurred. 
Comparing Figures [3] and [5] with the bottom part of Sudol & Harvey's (2005) Figure 1, 
the profiles of 5^^, B^^ and B^^ in Figures js] and [5] have clear, stepwise changes with no 
sign of an emission artifact. This serves as evidence that, while emission transients may have 
aflPected some of the pixels, the calculations of the integrated quantities plotted in Figures |3}(6] 
were not significantly compromised by artifacts. 

Figure [7| shows the temporal behavior of the average tilt angle for each fiare, the integral 
of 7 = tan~-^([S^ + B'^Y^'^ /\Br\) over the boxes around the neutral lines shown in Figure [l] 
divided by the area of the box in each case. Also shown is the field-strength-weighted average 
7 of 7, of the form 



^= B-fdA / B dA. (2) 

The changes in the tilt angle for the six fiares, and other magnetic changes discussed below, 
are summarized in Table [2| In all cases the tilt angle increased on average (the field became 
more horizontal), except, arguably, for the 2011 September 7 XI. 8 fiare. Figures [3] and [4] 
show that B^^ and B^^ both increased during the 2011 September 7 XI. 8 fiare, making the 
increase in the tilt profile less strong for this fiare than for the other fiares, whose vertical 
field changes are not so significant. The increasing pre-fiare trend of B^^ and the decreasing 
pre-fiare trend of B^^ give the tilt profile a steep decline before the fiare, that shows the 
lack of real significance of the tilt increase during this fiare. 

The temporal behavior of the average angle a = tan~-^(|5|| |/|S^|) between the horizontal 
field vector and the direction normal to the neutral line (approximated by the orientations 
of the short edges of the boxes around the neutral lines in Figure [T]) are shown in Figure [s} 
Also shown is the field-strength- weighted angular displacement 5, calculated in a manner 
analogous to Equation [2] In the plots in Figure [8| zero angular displacement would be 
consistent with an unsheared arcade whereas 90° would represent field aligned with the 
neutral line. The azimuthal angular displacement has been identified by some authors (e.g., 
Aulanier et al. 2012) with magnetic shear. This is because an unstressed magnetic field is 
expected to cross the main neutral line at an angle approximately normal to the neutral line. 
During the four largest fiares the horizontal neutral-line fields abruptly became less aligned 
with the normal direction, i.e., more aligned with the neutral line. The angular changes 
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Fig. 7. — Shown here are the average tilt angles (solid lines) and the field-weighted average 
tilt angles (dashed lines) near each neutral line plotted against time. These angles are 
averaged over the areas of integration indicated by the black rectangles in Figure [T] The 
vertical lines represent the GOES fiare start, peak and end times. 
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Fig. 8. — Shown here are the average azimuthal angular displacements from the normal to 
the neutral lines (solid lines) and the field- weighted average azimuthal angular displacements 
(dotted lines), at each neutral line plotted against time. These quantities are averaged over 
the areas of integration indicated by the black rectangles in Figure [TJ The vertical lines 
represent the GOES fiare start, peak and end times. 



16 



range from a few to nearly 10°. The 2011 March 9 XI. 5 flare flelds show some evidence of 
the same pattern but there is no such signature above the background evolution for the 2011 
February 13 M6.6 flare. 

Based on the evidence of Figures [7| and [8| the flelds show a strong tendency to become 
more tilted and more aligned with the neutral line as a result of flares. This is because the 
main effect of the flare on the photospheric neutral-line fleld is to strengthen the horizontal 
component parallel to the neutral line. 



3. Comparison with potential fields 

Although potential flelds do not accurately model the highly fluid-dominated magnetic 
fleld of the solar photosphere, they do represent Maxwell-stress- free states that are useful for 
diagnosing the stresses in observed non-potential flelds. We compute the shear of the mea- 
sured photospheric fleld with respect to a reference potential fleld as did Wang et al. (1992), 
Gosain and Venkatakrishnan (2012), Wang et al. (2012) among other authors. The full and 
horizontal magnetic shear angles are deflned by. 



S = cos-^[(B.BP)/(SSP)], (3) 
= cos-^[(Bh-B^)/(5h5,P)], (4) 

where B^ is the unique potential fleld such that = B^. In equations ^ and Q, 5^ = 
[{B?f + (5^)2 + (5P)2]i/2, Bl = {B%-Bl) and Bl = + (5P)2]1/2. The temporal 

behavior of the average full and horizontal shear angles for each flare is plotted in Figures [9] 
and [T0| derived by averaging S and S\, over the boxes around the neutral lines shown in 
Figure [T} Also shown are the fleld-strength-weighted full and horizontal shears, calculated 
in a manner analogous to Equation Q. The average of the magnetic shear = BS is also 
plotted in Figure [9j showing that this quantity is dominated by the behavior of S during 



every flare. From Figures M and 10 it is clear that the full shear angles are smaller than 



the horizontal shear angles, but the changes in are less signiflcant than the changes in S. 



The plots in Figures [9] and [10] do not follow the dominant patterns of the azimuthal angular 
displacements in Figure [8| There is some evidence that the full shear increased abruptly 
during the 2011 February 15 X2.2 and 2012 March 7 X5.4 flares. Three of the other flares 
did not show a signiflcant, abrupt change in shear. The 2011 September 7 XI. 8 flare fleld 
showed a marked permanent decrease in shear, particularly horizontal shear. However, this 
decrease began signiflcantly before the flare start time and ended after the flare end time. 
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Fig. 9. — Shown here are the average shear angles (sohd hnes) and the field-weighted average 
shear angles (dashed lines) near each neutral line plotted against time. The average of the 
magnetic shear is also plotted (dotted lines) in units of G°, divided by 1.5 x 10^ for the 2011 
February 15 and March 9 fiares and 10^ for the other fiares. These quantities are averaged 
over the areas of integration indicated by the black rectangles in Figure [ij The vertical lines 
represent the GOES fiare start, peak and end times. 
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Fig. 10. — Shown here are the average horizontal shear angles (solid lines) and the field- 
weighted average horizontal shear angles (dashed lines) near each neutral line plotted against 
time. These quantities are averaged over the areas of integration indicated by the black 
rectangles in Figure [T] The vertical lines represent the GOES fiare start, peak and end 
times. 
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Fig. 11. — Shown here are the average vertical dip angles near each neutral line plotted 
against time. These angles are averaged over the areas of integration indicated by the black 
rectangles in Figure [Tj The vertical lines represent the GOES flare start, peak and end times. 
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making it difficult to attribute the decrease in shear to the flare alone. In fact at the height 
of the flare the steep decline of the shear may have been interrupted somewhat, but there 
seems to have been no signiflcant increase of shear. The stepwise decrease in shear around 
the time of this flare appears to contradict the increasing azimuthal angular displacement 
show in Figure |8j As seen in Figure [2| the fleld changes resulted in enhanced vertical flelds 
on both sides of the neutral line that were distributed asymmetrically on both sides of the 
neutral line: on the north (positive) side the enhancements occurred at the east end of the 
neutral line while they occurred mostly at the west end on the south side. As a result, 
the potential fleld associated with the pre-flare measurements crosses the neutral line at 
a greater angle than the potential fleld based on measurements taken after the main fleld 
changes took place. Thus, even though the azimuthal displacement increased during this 
flare, the difference between the measured and potential horizontal flelds actually decreased 
on average and so the shear angle decreased. In summary, the only signiflcant and permanent 
shear changes directly attributable to the flares are the shear increases occurring during the 
2011 February 15 X2.2 and 2012 March 7 X5.4 flares. As Table [2] shows, the changes in 
azimuthal displacement show a more consistent pattern than the shear changes. 

We next compare the difference A7 = 7 — 7P between the observed tilt angle 7 and the 
potential fleld tilt angle 7^ = tan"^([(5^)^ + We refer to these parameters as 
the average and fleld-strength- averaged dip angles. The average and fleld-strength- weighted 



average of A7 are plotted against time for each flare in Figure [TT] A negative dip angle 
represents a fleld B more vertical on average than the reference fleld B^. These parameters 
clearly show uniform behavior over the whole data set. In each case the dip angle abruptly 
increased by several degrees during the flare, from a negative value to a value close to zero. 
During the 2011 September 6 X2.1 flare the dip angle abruptly changed from about —4° 
to about 2° before relaxing close to 0° over the next flve hours. Before the 2012 March 7 
X5.4 flare the dip angle was between —15° and —20° and it abruptly rose to about —5° 
during the flare. Over the flve-hour period before the 2011 February 13 M6.6 flare the 
dip angle steadily decreased from about —1° to about —6°, before snapping back to about 
—2° during the flare. In all cases, therefore, the pre-flare neutral-line fleld was signiflcantly 
less tilted (more vertical) than the corresponding potential fleld and, during each flare, this 
stress on the fleld was suddenly relaxed and the observed tilt angle more closely matched 
the tilt angle of the potential fleld. This pattern is consistent with the results of Gosain and 
Venkatakrishnan (2010) based on Hinode/SOT vector data for the 2006 December 13 X3.4 
flare. 

As Table [2] shows, these results for the dip angles are more uniform than the tilt angle 
proflles in Figure [7| whose left panels show more modest signs of change. The September 7 
XI. 8 flare tilt angle proflle is complicated by the increase of vertical flux shown in Figure [4} 



- 21 - 



The dip angle profile of this flare shows a much clearer stepwise change. This is because the 
sudden increase of vertical flux during the flare causes the potential fleld to become more 
vertical, with tilt angle closer to the tilt angle of the observed fleld. The abrupt relaxation 
of the fleld tilts towards the potential-fleld values is consistent with the interpretation of 
Hudson, Fisher and Welsch (2008) the coronal fleld reconnecting, relaxing and contracting 
into a simpler and more tilted conflguration. 



4. The electric current 



Figure 12 shows the evolution of the integrated vertical current flux, J^^, near the 
neutral line for each flare. The current proflles in the flgure are rather diverse but there 
are similarities between the electric current proflles of the two flares occurring in AR 11158, 
the 2011 February 13 and 15 flares, and also between those occurring in AR 11283, the 

2011 September 6 and 7 flares. During both the 2011 September 6 and 7 flares the negative 
current increased abruptly and permanently, creating a more balanced current in one case 
and a more imbalanced current in the other. These changes appear to be related to the 
increases in positive vertical magnetic flux during both flares, shown in Figure [4| Also the 
positive current decreased in a stepwise manner during the September 6 flare. The negative 
current decreased abruptly during the 2012 March 7 X5.4 flare. These last two changes were 
more signiflcant than any changes in the vertical magnetic flux (Figure [4]) . The vertical 
current reached local maxima during the 2011 February 13 and 15 flares but there was only 
an abrupt, permanent change in the negative current during the February 15 flare. These 
current maxima were not accompanied by corresponding peaks in the magnetic flux proflles 
in Figure [4] and so may represent magnetic stresses that were built up before the flares and 
released during them. A similar interpretation may apply to the current reduction during the 

2012 March 7 X5.4 flare. On the other hand, during the flares when the electric current was 
observed to increase abruptly near the neutral line of AR 11283, the 2011 September 6 and 7 



flares, the full and horizontal magnetic shear angles (Figures |9] and 10) either did not change 
permanently (September 6) or actually decreased (September 7). The interpretation of the 
magnetic changes during these flares is complicated by the abrupt and permanent vertical 
flux increases seen in Figure [4] which changed the reference potential fleld signiflcantly as 
discussed in Section [3| This is likely related to the fact that the vertical magnetic flux 
increased near the neutral line during both flares. Recall from Section |3] that the potential 
fleld associated with the September 7 pre-flare measurements crossed the neutral line at 
a greater angle than the potential fleld based on measurements taken after the main fleld 
changes took place, accounting for the reduction in magnetic shear there. The increase in 
current during the September 6 and 7 flares may therefore be mostly due to the increase in 
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Fig. 12. — The integrated vertical electric current flux near each neutral line, J^^ is plot- 
ted here for each flare as a function of time. The red/blue solid/dashed lines represent 
positive/negative current. The vertical lines represent the GOES flare start, peak and end 
times. The areas of integration for the neutral-line calculations are indicated by the black 
rectangles in Figure [TJ 
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field strength near the neutral line of AR 11283 during each flare. Finally, the March 9 flare 
appears not to have a significant electric current signature. It is clear that the flare-related 
electric current changes, like the shear changes, do not follow a general pattern. 



5. The Lorentz force changes 

We use the results of Fisher et al. (2012) to estimate the changes in the Lorentz force 
vector acting on the volume below the photosphere as a result of each flare. Assuming that 
the photospheric vector field is observed over a photospheric area at two times, t = 
before the field changes begin, and t = St after the main field changes have occurred, the 
corresponding changes in the Lorentz force vector components between these two times are 
given by Equations (17) and (18) of Fisher et al. (2012): 



and 




where at a fixed location in the photosphere 

SBl = BliSt)- 
5Bl = Bl{5t)- 

5{BrB^) = Br{St)B 



- SBl) dA, (5) 

(5,Bh) d^, (6) 

BliO) , (7) 

Bm , (8) 

(5t)-S,(0)Bh(0) . (9) 



The Lorentz force acting on the atmosphere above the photosphere is equal and opposite to 
the force acting on the volume at and below the photosphere (Fisher et al. 2012). 



Figure [T3| shows for each flare a spatial map of the Lorentz force changes and F^, 
derived by evaluating the integrals in Equations ([5|) and ^ pixel by pixel. The sums of 
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Fig. 13. — The Lorentz force vector changes during each flare. The vertical component, SF^ 
is indicated by the color scale and the horizontal components by the arrows with saturation 
values 10^ dynes/cm^ for the color scale and 2.5 x 10^ dynes/cm^ for the arrows. Red/blue 
coloring represents positive/negative (upward/downward) Lorentz force change. The black 
rectangles mark the regions of major field change near the neutral lines. The solid and dotted 
contours indicate strong {\Br\ > 1000 G) and quite strong {\Br\ > 100 G) fields, respectively. 
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the distributions shown in Figure [13] over the entire photospheric area gives the estimate for 
the total Lorentz force vector described by Fisher et al. (2012). Petrie (2012) argued that 
this calculation also gives a useful estimate of the spatial distributions of the Lorentz force 
vector components across the region for well-resolved changes in major magnetic structures. 
Alvar ado- Gomez et al. (2012) gave an analogous argument for the corresponding volume 
integral within the solar interior. 



Figure [13] shows the spatial distributions of the changes in the Lorentz force components 
during the flares. As we saw in Figures [4] and [5} the horizontal fleld changes (^B^ near the 
main neutral line increased the horizontal fleld strength, SB^ > 0, and were signiflcantly 
greater than the vertical fleld changes 5B'^. Equation [5] therefore leads us to expect the 
vertical Lorentz force changes to have been predominantly downward near the neutral lines. 



The plots in Figure 13 conflrm that near the main neutral lines the Lorentz forces clearly 
acted downwards into the photosphere during each flare. This behavior was anticipated to 
occur near neutral lines of flaring active regions by Hudson, Fisher and Welsch (2008) and 
Fisher et al. (2012), and has been found in past estimates of Lorentz force changes by Wang 
and Liu (2010), Petrie and Sudol (2010) and Petrie (2012). 

We now analyze the horizontal Lorentz force changes in a similar fashion. According 
to Equation ([6]), wherever the vertical fleld does not change signiflcantly compared to the 
horizontal changes and is positive/negative, the horizontal Lorentz force changes 5Fh should 
be parallel/anti-parallel to the horizontal fleld changes 5Bh. We already know from Figure [2] 
that on both sides of the neutral line 5Bh pointed approximately parallel to the neutral line 



in all cases except the 2011 March 9 flare whose pattern is not so clear. Figure 13 shows 
that, during the other flve flares, the Lorentz force acted in opposite directions on each side 
of the neutral line, with the changes on the positive side pointing parallel to 5Bh and those 
on the negative side anti-parallel to 5Bh as expected. 



Figure 14 shows the Lorentz force changes SF^ in the vertical direction as functions 



of time. The red and blue curves in Figure 14 represent positive (upward) and negative 



(downward) force changes, respectively. They were derived by taking running differences 
between consecutive image pairs and were integrated over the areas Anl represented by the 
black rectangles shown in Figures [T} [2] and [T3| The running difference plots in Figure 14 



show sharp, spiked signatures of the abrupt downward Lorentz force changes well above the 
noise levels during all six flares. It is clear from Equations ^ and ^ that the Lorentz force 
changes should resemble the time-derivatives of the fleld changes. Step-like fleld changes of 
brief but flnite duration, such as those seen in Figure [T3| produce the brief spikes of Lorentz 
force change in Figure [14} The fact that the dominant downward spikes in these plots are 
not matched by upward changes implies that the force changes were permanent. The black 
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Fig. 14. — The Lorentz force vector components in the vertical direction, SF^^^ are plotted 
as functions of time. The positive and negative running-difference changes are represented 
by the red dotted and blue dashed lines, respectively. The fixed-difference changes with 
respect to the first image are represented by the solid black lines. The areas of integration 
are indicated by the black rectangles in Figure [Tj The red/blue dotted/dashed lines represent 
positive/negative force changes. The vertical lines represent the GOES fiare start, peak and 
end times. 
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Fig. 15. — The Lorentz force vector components in the paraUel direction, SF^^^ are plotted 
as functions of time. The red sohd and blue dashed lines represent the force changes acting 
on the positive and negative sides of the neutral lines, respectively. The areas of integration 
are indicated by the black rectangles in Figure [!} Positive/negative force changes act in 
the approximately westward/eastward directions along the long edges of these rectangles, 
respectively. The vertical lines represent the GOES flare start, peak and end times. 
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Fig. 16. — The Lorentz force vector components in the perpendicular direction, SF^^^ are 
plotted as functions of time. The red solid and blue dashed lines represent the force changes 
acting on the positive and negative sides of the neutral lines, respectively. The areas of 
integration are indicated by the black rectangles in Figure [l} Positive/negative force changes 
act in the approximately northward/southward directions along the short edges of these 
rectangles, respectively. The vertical lines represent the GOES flare start, peak and end 
times. 
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curves in Figure [14] show force changes derived by fixed differences with respect to the first 
image. These curves confirm that aU of the fiare-related force changes were permanent. They 
also show evidence of a steady net upward Lorentz force before three of the fiares, the 2011 
February 15 X2.2, September 6 X2.1 and September 7 XI. 8 fiares, perhaps associated with 
pre-fiare infiation of the fields near the neutral lines. Gosain (2012) described a pre-fiare rise 
phase in a study of AIA observations of the 2011 February 15 X2.2 fiare. 

The sizes of these force changes, ranging from about 9 x 10^^ dynes to about 6 x 
10^^ dynes, are comparable to those found in the previous estimates of fiare-related Lorentz 
force changes by Wang and Liu (2010) and Petrie and Sudol (2010). For the 2002 July 
26 M8.7 fiare Wang and Liu (2010) found a downward force change of 1.6 x 10^^ dynes. 
Petrie and Sudol (2010) found a range of longitudinal force change estimates up to about 
2 X 10^^ dynes. Petrie and Sudol's estimates are likely to have been underestimates because 
they included only information on the longitudinal field component. 



Figures 



15 



and 



tal directions paral 



16 



show the Lorentz force vector changes, 6F^^ and 5F^^, in the horizon- 
el and perpendicular to the neutral line as functions of time. The forces 
are integrated over the areas ^nl represented by the black rectangles shown in Figures [T} [2] 
and 13, The parallel directions are the directions of the long edges of the rectangles, point- 
ing approximately west. The perpendicular directions are the direction of the short edges 
of the rectangles, pointing approximately north. Figures [15] shows significant and abrupt 
force changes at the time of every fiare. The red and blue curves describe the parallel force 
changes occurring on the positive and negative sides of the neutral line, respectively. In each 
case the force changes are directed in opposite directions on the two sides of the neutral line. 
The corresponding perpendicular force changes, shown in Figure [16} were smaller and less 
significant, and were generally not permanent. 



Table [3] summarizes the directions of the parallel and perpendicular Lorentz force 
changes, SF^^ and associated with the six fiares. In all cases where permanent 

changes of Lorentz force were detected on the northern/southern side of the neutral line, 
the perpendicular force change (columns 7 and 8 of Table js]) was directed south/north, i.e., 
towards the neutral line. Both sides of the neutral line were observed to undergo significant, 
permanent perpendicular force change during only one fiare, the 2011 March 9 XI. 5 fiare. 
There were no significant, permanent changes on either side of the neutral line during the 
2012 March 7 X5.4 fiare, and one side underwent permanent change during the four other 
fiares. The parallel force changes (columns 5 and 6 of Table js]) were more significant, and 
were detected in both polarities during all fiares except the 2011 March 9 XI. 5 fiare, when 
the north side of the neutral line underwent significant, permanent change while the south 
side didn't. In each case when both sides of the neutral line experienced permanent paral- 
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lei Lorentz force changes, these occurred in opposite directions on each side of the neutral 
line. In columns 5 and 6 of Table |3] the directions of the positive polarities are in bold 
typeface. These directions agree perfectly with the axial field directions (column 4). This 
implies that the positive polarity footpoints of the fields near the neutral line were pulled 
towards the corresponding negative footpoints. Meanwhile, the negative footpoints were 
being pulled towards the positive footpoints during five of the six fiares, and underwent no 
permanent change during the remaining fiare. Ignoring vertical movements, a convergence 
of footpoints along the direction parallel to the NL with no movement in the perpendicular 
direction would produce a decrease of shear, whereas a convergence of footpoints along the 
direction perpendicular to the NL with no movement in the parallel direction would produce 
an increase of shear. These interpretations ignore the effects of vertical movements such as 
loop collapse. The parallel force changes alone would therefore have relaxed the shear of 
the neutral line fields, whereas the perpendicular forces acted towards the neutral line and 
would have tended to increase the shear. 

The perpendicular force changes were not generally significant compared to the back- 
ground evolution. Furthermore, comparing the sizes of the force changes parallel and per- 
pendicular to the neutral line, the ratio SF^^/SF^^ is less than 0.1 for the large fiares and 
less than 0.5 for the small fiares. For comparison, the ratio Bf^/Bj^^ ranges from about 
0.3 to about 0.5 (Figures [5] and [6]). The horizontal force changes therefore can't explain the 
changes of azimuthal diflFerence and shear angles seen in Figures [8| [9] and [T0| They likely 
represent the contraction of the collapsing, line-tied loop fields, trying to pull the two sides 
of the neutral line towards each other. 



6. Conclusion 

In 12 hour time series of 12-minute SDO/HMI vector field observations covering six 
major fiares we have found consistent patterns in the field and Lorentz force vector changes 
near the main magnetic neutral lines of the fiaring regions. We summarize the main results 
before drawing conclusions from them. 

1. Near the main magnetic neutral lines, the field vectors became stronger and more hor- 
izontal during all six fiares. This was almost entirely due to an increase in strength of 
the horizontal field components parallel to the neutral line in each case. The horizontal 
perpendicular and vertical components did not show comparably significant and per- 
manent changes during the fiares. Even during the 2011 September 7 fiare, when the 
vertical field did change significantly and permanently, the change in the horizontal 
parallel field was more significant. 



-31 - 



2. The magnetic tilt angle increased significantly during all flares with the arguable ex- 
ception of the 2011 September 7 flare. The neutral line flelds were less tilted than 
their corresponding potential flelds before all six flares and relaxed abruptly and per- 
manently closer to potential fleld tilt angles during every flare. 

3. The horizontal flelds became signiflcantly and permanently more aligned with the 
neutral line during the four largest flares. The (full and horizontal) shear angle with 
respect to the reference potential fleld increased signiflcantly during only two of the 
flares. 

4. The electric current near the neutral line showed a marked change of behavior around 
the times of most of the flares but there was no consistent pattern in the character of 
this change. 

5. The vertical Lorentz force had a large, abrupt, permanent downward change during 
each of the flares. The sizes of these force changes, ranging from about 9 x 10^^ dynes 
to about 6 X 10^^ dynes, are comparable to those found in the previous estimates of 
flare-related Lorentz force changes. 

6. The horizontal Lorentz force component parallel to the neutral line showed signiflcant, 
permanent changes acting in opposite directions on each side of neutral line during 
each of the six flares, comparable in size to the vertical force changes. In all cases the 
shearing forces were consistent with a decrease of shear near the neutral line, whereas 
the fleld itself became more aligned with the neutral line as a result of the increase 
in the horizontal fleld strength. The horizontal force changes perpendicular to the 
neutral line tended to act towards the neutral line but these were not as signiflcant as 
the parallel changes. 

The abrupt and permanent increase in magnetic fleld strength near the neutral line 
during a flare, accompanied by corresponding increase in magnetic tilt and downward Lorentz 
force change, is most easily explained by the flux near the neutral line being compressed 
from above. Because the changes in the horizontal fleld were generally not accompanied by 
corresponding changes in the vertical fleld, the resulting changes of fleld tilt cannot have 
been the result of simple rotation of the magnetic vector towards the neutral line but must 
have come from the collapse of nearly horizontal magnetic fleld towards the photospheric 
magnetic neutral line from the surrounding volume. In each case the pre-flare fleld was less 
tilted than the reference potential fleld and abruptly collapsed during the flare to a tilt angle 
much closer to that of the potential fleld. This implies that, before each flare, the fleld was 
supported in a more vertical conflguration and relaxed closer to a potential-fleld tilt angle 
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during each flare. The large downward Lorentz force changes observed during all six flares 
are consistent with the interpretation of sheared loop collapse. 

We flnd that the tilt angles of the observed flelds with respect to the potential flelds 
has a stronger increasing pattern than the tilt with respect to the vertical direction, whereas 
the behavior of the azimuthal angles of the observed flelds is more consistent with respect to 
the neutral line than with respect to the potential flelds. The observed fleld change during 
each flare was dominated by an abrupt, permanent and signiflcant increase in the horizontal 
component parallel to the neutral line, and horizontal fleld became signiflcant ly more aligned 
with the neutral line during the four largest flares in the study. These results imply that 
the collapsing fleld was almost aligned with the neutral line in each case, and more aligned 
on average with the neutral line than the pre-flare neutral line fleld. On the other hand, 
the average shear angle increased during two flares and decreased during one, demonstrating 
that the flare process does not generally relieve magnetic stresses associated with magnetic 
shear. The lack of a general pattern in the vertical electric current changes also reflects this. 
The change in parallel horizontal flelds and increased alignment between neutral line flelds 
to neutral lines appear simply to be consequences of the relaxation of stresses associated 
with the non-potential tilt angles. 

Wang (2006) found an unshearing movement parallel to the neutral line in flare-related 
longitudinal magnetic fleld changes in all flve events that he studied, consistent with overall 
release of shear. In our results, the un-shearing pattern of the generally dominant parallel 
horizontal Lorentz force changes would by themselves have reduced the shear of the fleld 
near the neutral lines during each flare, but the magnetic shear only decreased signiflcant ly 
during one of the six flares. These un-shearing patterns can be interpreted as a signature of 
the neutral line flelds contracting during the flare, pulling the two sides of the neutral line 
towards each other. It seems that the horizontal force changes did not prevail in un-shearing 
the neutral line flelds because the line-tying effect of the dense photospheric plasma makes it 
much more difflcult to move photospheric footpoints laterally than to change their tilt angle. 
This may explain the differing physical effects of the horizontal and vertical Lorentz forces: 
vertical, tilt-related magnetic stresses appear to be much more easily relieved by flares than 
horizontal, shear-related stresses. 

Nonlinear force-free fleld extrapolations have presented a consistent picture of a sheared 
structure collapsing towards the neutral line, leaving a void above that is flUed by more 
relaxed fleld (Jing et al. 2008, Sun et al. 2012, Liu et al. 2012). This modeling work suggests 
that if magnetic shear increases at low altitudes during flares it may decrease above a certain 
height, producing a net decrease of shear in the system as a whole. Our results suggest that 
the main flare-related fleld changes are caused by the release of magnetic stresses associated 
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with non-potential tilt/dip angles, and that the consequent strengthening of the parallel 
horizontal field component may or may not increase photospheric magnetic shear. The 
scenario of sheared loop collapse developed above would lead to a reduction of shear higher 
in the atmosphere in agreement with the modeling results. 

We have only studied the effects of major fiares on fields near neutral lines in this 
paper. While clear patterns have emerged from this sample, it will be instructive to examine 
a large sample of high-cadence vector-field measurements from HMI and the National Solar 
Observatory's Synoptic Optical Long-term Investigations of the Sun (SOLIS) telescope to 
discover if these patterns are generally dominant, and whether an analogous pattern can be 
found in the fiare-related magnetic changes observed in sunspots. 

I thank the referee for a stimulating review that helped to develop the project. I thank 
Sanjay Gosain and Alexei Pevtsov for discussions. SDO is a mission for NASA's Living With 
a Star program. This work was supported by NSF Award No. 106205 to the National Solar 
Observatory. 
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